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By means of space-group symmetry arguments, we argue that the electronic pairing in iron-based high temper- 
ature superconductors shows a structure which is a linear combination of planar s-wave and d-wave symmetry 
channels, both preserving the 3-dimensional A\g irreducible representation of the corresponding crystal point- 
group. We demonstrate that the s- and d-wave channels are determined by the parity under reflection of the 
electronic orbitals through the iron planes, and by improper rotations around the iron sites. We provide evidence 
of these general properties by performing accurate quantum Monte Carlo ab-initio calculations of the pairing 
function, for a FeSe lattice with tetragonal experimental geometry at ambient pressure. We find that this picture 
survives even in the FeSe under pressure and at low temperatures, when the tetragonal point-group symmetry 
is slightly broken. Our theory can rationalize and explain a series of contradictory experimental findings, such 
as the observation of twofold symmetry in the FeSe superconducting phase, the anomalous drop of Tc with 
Co-impurity in LaFeAsO(i_a;)Fi:, the s-to-d-wave gap transition in BaFe2As2 under K doping, and the nodes 
appearing in the LiFeAs superconducting gap upon P isovalent substitution. 
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The pairing symmetry of the superconducting state in iron- 
based superconductors (IBS) has been one of the most de- 
bated subjects since their first discovery, in both theory and 
experiments. ifH-llt] Its determination is paiticulaiiy challeng- 
ing in the IBS for their complex electronic structure, with a 
strong multiband character, and a Fermi suiface constituted 
by many sheets, which can vary with doping and chemical 
composition. Moreover, the IBS families are usually compen- 
sated metals, and thus both electrons and holes are involved 
in the paired state. Weak-coupling RPA approaches, cou- 
pled to multiband BCS theory, yield a pairing function with 
global s-wave {Aig) symmetry, but with electron and holes 
pockets having opposite sign. This scenario, dubbed s*, was 
first proposed by Mazin et al.\A\ complemented later with its 
variants, called "extended s^",|5] in the case that accidental 
nodes appear on a single Fermi sheet without breaking the 
full rotational symmetry. The latest generalizations include 
also "weak" nodal lines which develop as closed loops on the 
3-dimensional (3D) Fermi surface. [6,-i8J 

A variety of experiments has been performed to probe the 
pairing symmetry of IBS, ranging from thermal conductivity 
and specific heat measurements to angle-resolved photoemis- 
sion spectroscopy. While there is no doubt on the spin singlet 
nature of the paiiing state as revealed by the Knight shift, 1 9] 
the presence of nodes and the total symmetry of the spatial 
part of the paiiing function are controversial. In fact, the ex- 
perimental outcome seems to lack universality, as it depends 
crucially on the "family" of tested compound, its doping, its 
isovalent substitutions, and its level of disorder For instance, 
the situation is paradoxical for the "111" family, where the 
LiFeAs material shows a fully gapped superconducting or- 



der parameter, while the isovalent substitution by phosphorus 
leads to thermodynamic properties compatible with a nodal 
paiiing function.! 10] For the "122" family, there is a recent 
claim, supported by independent experimental probes, that 
the pairing in the BaFe2As2 undergoes an s-Xo-d symmetry 
change by doping with potassium,| 1 1] the fully substituted 
KFe2As2 being identified as a d-wave superconductor [, 1 7i In 
the "1111" family, the dependence of the critical temperature 
Tc upon disorder has been studied in the LaFeAsO(i_j.)Fa; 
with X = 0.11. It has been found that Co-doping induced dis- 
order makes Tc to fall much more slowly than Mn-doping.lll3l 
This is not compatible with any sign changing order param- 
eter (either or d-wave). Finally, specific heat measure- 
ments performed on the FeSe revealed a highly anisotropic 
order parameter in the "11" family, with a twofold sym- 
metry directly observed by scanning tunneling microscopy 
(STM).|15] 

In this Letter, we study the IBS gap structure from a dif- 
ferent perspective, namely by looking for its universal fea- 
tures based on symmetry constraints induced by the 3D space- 
group transformations of the crystals. We prove that the IBS 
pairing function is a linear combination of terms having pla- 
nar s- and d-wave symmetries, both fulfilling the full 3D Aig 
representation. These properties are then verified in the FeSe 
by performing state-of-the-art quantum Monte Carlo (QMC) 
calculations from first-principles. Our theory provides a gen- 
eral framework to account for the contradictory experimental 
outcomes. 

The undistorted structure of FeSe, together with the other 
parent compounds of the "111" and "1111" families, belongs 
to the PA/nmm space-group, which is nonsymmorphic due 



2 



to a glide plane parallel to the iron layers. Indeed, its related 
point-group classes contain certain operations with the 
nonprimitive lattice translation r = (1/2, 1/2, 0) - expressed 
in crystal units of the 2-Fe unit cell - which brings a Fe site 
into the other in the cell. This is due to the Se sites alternating 
above and below the Fe layer In particular, the C4 class of 
pure rotations by 7r/2 about the principal axis is not compati- 
ble with the 3D lattice, unless associated with the translation 
T. On the other hand, the iS'4 class, defined as a rotation by 
tt/2 about the principal axis combined with the reflection ah 
through the Fe plane (S4 = C^ah), is a a pure point-group op- 
eration allowed by the PA/nmm space-group, as well as the 
reflection through the plane perpendicular to the Fe (xy) plane 
and bisecting the secondary (x and y) axes (the (T^ class). 

This has important consequences on the symmetries of the 
subgroup associated with the 2-dimensional (2D) Fe square 
sublattice, and thus on the low-energy features of the system, 
as for instance, the superconducting gap structure. Indeed, the 
standard C^i, group of the 2D square lattice, generated by the 
elementary (Jd and C4 point-group operations, does not ful- 
fill all symmetries of the PA/nmm space-group, whereas the 
group generated by <Td and ^4 does. Moreover, proper (C4) 
and improper (S4) rotations differ from the physical point of 
view, because a Hamiltonian or a pairing function fulfilling 
the point-group operations generated by ^4 and (Td will trans- 
form in a very different way with respect to those based on the 
standard C41, group, as it will be shown in the following. 

From ab-initio electronic structure DFT and DMFT calcu- 
lations of iBs. fiemsli it is well established that the bands 
crossing the Fermi level have a strong d-orbital character 
Indeed, the minimal tight-binding Hamiltonian which repro- 
duces successfully the Fermi surface contains all the five 
atomic (or Wannier) d-orbitals centered on the Fe sites. 11 1911 
namely \d^^), jdyz), \dxy), \d^2_y2), and \d3^2_^.2), ex- 
pressed with the crystal axes oriented according to the Fe 
square lattice directions. A given atomic orbital (iR,iy.(j(r) 
defining the 2D multiorbital model can be represented by 
(r|(iR,;y.(j) = ('"|cr (j|0}, where r lives in the 3D space, 
R is the center of the atomic orbital on the square lattice, h' 
is the orbital symmetry, and <t = ±-i is the electron spin. 
A spatial point-group operation acts on the orbital cre- 
ation operator by changing its position R i?,,R, its phase 
s = 1 ^ s'{ri, v) — ±1 and its orbital label v — >■ v) 
according to the general rule: 



For rj — {ad, 54, C4} and v given, the corresponding s' and v' 
are reported in Tab. I] The transformation rules are different, 
whether the group is generated by ^4 or C4. 

Let us take into account the 2D tight-binding Hamiltonian 
Ha which defines the low-energy model for the IBS. It reads 
as 
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TABLE I: Action of the symmetry operations r; — {aa, 54, C4} on 
the five iron atomic d-orbitals. The entries represent the pair s' — i: 
and v' = dxz,dyz,dxy,d^2_y2,d^^2_^2 as a function of the or- 
bital index u (first column) and the symmetry operation rj considered. 
Note the different s' between proper C4 and improper S4 rotations 
in the first two orbitals (odd with respect to 0"^). 



where the vectors {k} belong to the Brillouin zone of the 
2D lattice, Ck,iy,CT (cj^j^o.) '■^^ Fourier transform of cr^i/^o- 
(c^j^^), the ^,1/ indexes run over the d-orbital labels, the 
hoppings fulfill the hermitian condition i(k)* ~ i(k)i,^, 
and t{—'k),j fj = i(k)* ^ holds for reality. Obviously, Hq must 
be invariant under the symmetry transformations of the group 
generated by ad and 6*4, as the compatible subgroup of the full 
PA/nmm space-group. If we apply the transformation rules 
in Eq. [T|to Hq, and require its invariance, we obtain a set of 
relations for the hopping matrix: 

(3) 

for the eight possible independent point-group operations, de- 
fined hy rj = ad, S4, and all their products. These transforma- 
tions have already been derived in Refs. 2o'and'2T. However, 
their implications for the pairing properties have so far been 
overlooked. 

The pairing operator F is assumed to be translationally in- 
variant (on the square lattice), so that it is convenient to write 
it in k-space as 



F = 



(4) 



(2) 



Singlet pairing, and no broken time-reversal symmetry imply 
that /(k)^,^ = /(-k)^,^, and /(k)^,^ = /(-k)*^, respec- 
tively. In the following, we derive the spatial-orbital proper- 
ties of the pairing function F(r, r') = (r, r'|F|0), by assum- 
ing that it belongs to the Aig irreducible representation of the 
point-group generated by ad and 5*4. This assumption will be 
verified later by ab-initio QMC calculations with a fully op- 
timized Jastrow projected BCS wave function. "Improper s- 
wave" is the name we assign to the symmetry of the resulting 
F, to distinguish it from the standard or "proper" one, i.e. the 
Aig irreducible representation of the isomorphic point-group 
where the S'4 class has been replaced by C4. 

Based on the general rules in Eq. [T]and the entries of Tab. |I] 
we can derive the symmetry constraints for the improper s- 
wave. As the characters of the Aig representation are - by def- 
inition - identically one for all classes within the point-group, 
the symmetry transformations of /(k)^_^ read as 

f{Rrjk)^,f, = s(ry, iy)s{f], M)/(k)i.'(^,i.),^'(^,^), (5) 
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for 77 = ad, 5*4, and all their independent products. The 
above constraints reduce the variational freedom for the signs 
of F(r, r'). Accidental nodes must not break the general sym- 
metries in Eq. |5] For the improper s-wave, we demonstrate 
that F has a peculiar, symmetry selected, sign structure. To 
show this, let us divide the iron d-orbitals into two groups, 
namely the "even" (e) and"odd" (o) ones, according to their 
parity under ah- The \dxz) and \dyz) orbitals are odd, while 
\dxy), \dx'^-y-i), and \d^z'^-r'^) are even. Thus, F can be ex- 
panded as 



F — F^^ F^Q ~\- Fq, 



Fo, 



(6) 



Foe(r,r') = 



with its components defined as follows: 

Fee(r,r') = Ek,{^,M}eeven/(k)<>,A.c'k.^,t(r)d_k,p,i(r'), 

I]k,{i/}eeven/(k)iy,pdk,i/,t(l')'^-k,Ai,4.(r')' 
{^i}eodd 

I]k,{i'}eodd/(k)^,^dk,,y,t(i")d-k,M,i(r')> 

{/i}£even 

Ek,{,.,p}Godd/(k)^,MC^k,,.,tWc^-k,p,i(r')- 

The above pairing functions are periodic in both r and r', 
with periodicity set by the Bravais lattice vectors. In the fol- 
lowing analysis, without loss of generality, we are going to pin 
r around an iron site [22], while we define F{kx,ky) as the 2D 
Fourier transform of F with respect to the planar (a;, y) com- 
ponents of r' — r. According to the transformation rules in 
Eq. |5] and in Tab. HI we obtain that Fee and Fqo have an s- 
wave planar symmetry, e.g. Fee{kx,ky) = Fee{-kx,ky) = 
Fee{kx,-ky) = Fee{ky,k^), whcrcas Feo and Foe, which 
couple unlike-parity components, have planar d^y symmetry. 



S-g- Feo{kx^ky^ — Feo(^ kx , ky 



Feo(^kx^ ky^ 



Feo{ky,kx). We note here that in a standard s-wave super- 
conductor obeying the C41, point-group transformations, also 
the Feo and Foe components are s-wave. Thus, our improper- 
s wave pairing function has to have different signs just from 
symmetry considerations. 

These predictions have been verified against accurate ab- 
initio QMC calculations based on the energy minimization 
of a Jastrow projected BCS wavefunction, performed for the 
FeSe on a 4 X 4 X 1 iron lattice, subject to periodic boundary 
conditions. We used the experimental FeSe geometry, deter- 
mined at ambient pressure (0 GPa) and temperature [l23ll and 
under a hydrostatic pressure of 4 GPa at low-temperature. [E4I1 
Fe and Se atoms have been replaced by pseudoatoms, con- 
taining only 16 and 6 valence electrons, respectively. ll25ll The 
full many-body wave function can be written in a compact 
form as *jBcs(R-ei) = exp[- J(Rei)] det[0(rj, rj)], where 
1 < ijj ^ N/2, with N the total number of electrons in the 
supercell, and Rei = {r|, . . . , r^^/^, • ■ • , many- 



body electron configuration. Dynamic correlations, including 
charge fluctuations, are taken into account by the Jastrow fac- 
tor J, while correlations built up in the superconducting state 
are described by the BCS pairing function (j), which reads as 



M 



(7) 




FIG. 1: Even-even (panel a), and odd-even (panel b) components of 
the pairing function 0'^°'"'" of the FeSe high-Tc superconductor at 
GPa obtained by QMC, by optimizing the energy of the variational 
wavefunction without assuming any global symmetry. The contour 



plots show 



■■(Re. 



, r) with Rc, 



■ set to be the iron lattice site 



at the center of the supercell, while r spans the plane defined by the 
4x4 lattice. Red (yellow) balls are iron (selenium) sites. Arbitrary 
units blue (red) intensity indicates negative (positive) regions with 
corresponding magnitude. 




FIG. 2: The same as in Fig.[T]but for a structure corresponding to a 
pressure of 4 GPa. 



The above expansion is on the lowest AI molecular orbitals 
(MO) V'f^'^j generated by density functional theory (DFT) 
calculations, and further optimized by QMC|25]. In all the 
calculations quantum effects on nuclei are neglected, within 
the Bom-Oppenheimer approximation. Hence superconduc- 
tivity is assumed to be non conventional, namely not coming 
from electron-phonon coupling. 

(/)(r, r') is an extension of F(r, r') in Eq.|6l as it is expanded 
over the full atomic basis set used to represent the MO's. The 
resulting (j>, restricted to its correlated part 0'=°'''" (obtained by 
summing the MOs in Eq. |2]over i with N/2 < i < M), 
is plotted in Figs. [T] and |2] for and 4 GPa, respectively. 
^corr jj^g been projected over its 4 possible components t/)™"' 
(Figs, ma) andEa)), ^^^^ and ^-f^^ (Figs.Eb) and 

|2b)), based on the parity with respect to a^ 1 26]. The fi- 
nal picture is remarkably consistent with the improper s-wave 
pairing function predicted by the 2D point-group symmetries 
and under the assumption of a 3D A\g irreducible representa- 
tion. Our QMC solution is A\g as a consequence of an unbi- 
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ased energy minimization, and it has not been imposed a pri- 
ori in the optimization. Thanks to the improper S4 generated 
point-group, the pairing function (p shows two planar symme- 
try channels, s-wave for the and (poo^'', dr^y-wawe for 
the and tp^"^^ components. We emphasize that (p does 
not possess a well defined symmetry under a proper rotation 
by 7r/2, just because such a rotation is not a symmetry of the 
crystal. Moreover, Fig. |2] shows that this picture holds even 
when the tetragonal point-group symmetry is slightly broken, 
and lowered to the orthorhombic point-group, as in the case 
of applied pressure and/or low temperature. 
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FIG. 3: Condensation energy in the 4x4 system (in Hartree), and 
pairing gap Abcs (in meV) at and 4 GPa, obtained from our QMC 
calculations, compared to the experimental superconducting temper- 
ature Tc in meV. 

To further check the reliability of 'I'jbcs to describe the su- 
perconducting features of FeSe, we estimated (via correlated 
sampling) the condensation energy €cond, defined as the dif- 
ference between the expectation value of H computed with 
the fully paired 5'jbcs and the one with the Jastrow pro- 
jected single determinant, obtained by summing the MOs with 
i < N/2 in Eq.|2] Another independent test has been done by 
fitting the variational parameters A; in Eq.|7]for i > N/2 with 
the BCS expression Abcs/ + ^{a - + ^scs)' 
where we took = ef^^ (the DFT energy of the correspond- 
ing MO), /i equal to the Fermi level, and A bcs the only fitting 
parameter Those quantities are much more sensitive to finite 
size errors than the pairing symmetry, set by the short-range 
behavior of (p, which is converged in the 4x4 supercell. The 
tcond and Abcs are shown in Fig.|3] together with the exper- 
imental critical temperature T^. It turns out that Abcs ~ Tc, 
which is a remarkable result given the tiny energy scale (of 
the order of 1 meV) involved in the gap opening. Both econd 
and Abcs follow quite well the behavior of the experimental 
Tc, which increases with pressure, in the pressure range con- 
sidered herel27]. Therefore, our pairing function describes 
reasonably well the experiments. Remarkably, we are close to 
make quantitative predictions on unconventional high temper- 
ature superconductors by means of a fully ab-initio method. 

We explore now the physical consequences of the improper 
s-wave from a general perspective, by assuming that also the 
other IBS families have a Aig superconducting gap, as the 
FeSe. First of all, the twofold symmetry of the pairing func- 



tion observed by STMfls'l in stoichiometric FeSe with impu- 
rities and magnetic field can be explained by assuming that 
the interference states probed by the STM tip are affected by 
the sign change induced by a planar superposition of the s- 
and d-wave channels, invariant under tt rotations. Note that 
there are two possible arrangements, oriented in a dxy fashion 
and tilted by 7r/2, in agreement with the different orientations 
measured by the STM zero bias conductance. 



From ab-initio DFT and DMFT calculations. 11 171 11811 the 



d,j.2_y2 and d^z^^r'^ orbitals are usually pseudogapped. This 
allows us to disregard them from the summation in Eq. |4] 
Thus, dxy is the only remaining orbital with even parity, and it 
can trigger the formation of the d-wave channel in the pairing 
function. Indeed, for absent or weak d^y spectral weight at the 
Fermi level, the two orbitals left in F are both odd under ah, 
and therefore they can only contribute to the s-wave channel. 

In the "111" family, DFT band structure calculations and de 
Haas-van Alphen experiments show that the main effect of P 
isovalent substitution is the formation of the third hole pocket 
- with dominant d^y character -, absent in the LiFeAs. l.28ll 
The enhanced d^y spectral weight in the P substituted com- 
pound and the subsequent development of an "even-odd" d- 
wave pairing gap can explain why the LiFeP pairing is nodal, 
while the LiFeAs shows a fully gapped behavior 

Within the improper s-wave picture, the properties of the 
"1111" family modified by impurity substitution as in Ref.ll3 
are consistent with a tendency of F to develop nodes by go- 
ing from Co to Mn doping through the opening of the d-wave 
channel. Indeed, Co impurities are electron donors, while Mn 
impurities dope the system with hole, thus enhancing the dxy 
spectral weight around F. 

The "122" family does not belong to the PA/nmra space- 
group. Nevertheless, by inspection of its lA/mmm space- 
group unit cell, one can see that also in that case the ^4 im- 
proper rotation is the valid symmetry on the iron layer, thus 
the improper s-wave theory applies also here. The s-to-d 
wave transition by doping the BaFe2As2 with K is one of the 
most striking experimental situations for us.[l IJ In the strong 
overdoped regime, the band structure certainly shows all three 
pockets at F with again an enhanced dxy spectral weight. If 
the d-wave channel dominates in the improper s-wave gap, 
then it is plausible that the experimental thermodynamic fea- 
tures are compatible with a pure d-wave behavior 

Finally, we would like to point out that the importance of 
the third hole pocket, with strong dxy character, has already 
been observed by Kuroki and co-workers in relation to the 
Tc dependence on the pnictogen height. [29J The improper s- 
wave symmetry explains the key role played by dxy, partic- 
ularly in connection with the appearance of gap nodes, the 
twofold anisotropy, and the stabilization of the d-wave chan- 
nel, by reconciling a series of experiments which look other- 
wise contradictory. It suggests also that the pnictogen or the 
chalcogen, rather than being spectators, take an active role in 
the electronic pairing, by bridging two next-nearest neighbor 
iron sites in a dxy fashion via the improper s-wave symmetry. 
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